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Abstract 

vo 

We outline the WIMP dark matter parameter space in the Constrained MSSM by performing a comprehensive statistical analysis 
"^h that compares with experimental data predicted superpartner masses and other collider observables as well as a cold dark matter 
(— ) abundance. We include uncertainties arising from theoretical approximations as well as from residual experimental errors on 
I relevant Standard Model parameters. 
O | We present high-probability regions for neutralino dark matter direct detection cross section, and we find that 1CT 10 pb < o^ 1 < 
* 10 -8 pb for direct WIMP detection (with details slightly dependent on the assumptions made). We highlight a complementarity 
^ between LHC and WIMP dark matter searches in exploring the CMSSM parameter space. We conclude that most of the 95% 
• • probability region for the cross section will be explored by future one-tonne detectors, that will therefore cover most of the currently 
> favoured region of parameter space. 
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1. Introduction 

Two of the most challenging questions facing particle 
physics today are the instability of the Higgs mass against 
radiative corrections (known as the "fine-tuning problem" ) 
and the nature of dark matter. Unlike the Standard Model 
(SM), weak scale softly broken supersymmetry (SUSY) 
provides solutions to both of them. Firstly, the fine-tuning 
problem is addressed via the cancellation of quadratic di- 
vergences in the radiative corrections to the Higgs mass. 
Secondly, assuming i?-parity, the lightest supersymmetric 
particle (LSP) is a leading weakly interactive massive par- 
ticle (WIMP) candidate for cold dark matter (CDM). De- 
spite these and other attractive features, without a refer- 
ence to grand unified theories (GUTs), low energy SUSY 
models suffer from the lack of predictivity due to a large 
number of free parameters (e.g. , over 120 in the Minimal 
Supersymmetric Standard Model (MSSM)), most of which 



arise from the SUSY breaking sector. 

The MSSM with one particularly popular choice of uni- 
versal boundary conditions at the grand unification scale is 
called the Constrained Minimal Supersymmetric Standard 
Model (CMSSM) [1]. The CMSSM is defined in terms of 
five free parameters: common scalar (mo), gaugino {m 1 / 2 ) 
and tri-linear (Aq) mass parameters (all specified at the 
GUT scale) plus the ratio of Higgs vacuum expectation val- 
ues tan/3 and sign(/x), where fi is the Higgs/higgsino mass 
parameter whose square is computed from the conditions 
of radiative electroweak symmetry breaking (EWSB). The 
economy of parameters in this scheme makes it a useful tool 
for exploring SUSY phenomenology 

Many studies have explored the phenomenology of the 
CMSSM or other SUSY models, mostly by evaluating 
the goodness-of-fit of points scanned using fixed grids in 
parameter space, see e.g. [2,3,4,5,6,7]. However, this ap- 
proach has several severe limitations. Firstly, the number 



of points required scales as k N , where N is the number 
of the model's parameters and k the number of points for 
each of them. Therefore this approach becomes highly in- 
efficient for exploring with sufficient resolution parameter 
spaces of even modest dimensionality, say N > 3. Sec- 
ondly, narrow "wedges" and similar features of parameter 
space can easily be missed by not setting a fine enough 
resolution (which, on the other hand, may be completely 
unnecessary outside such special regions). Thirdly, extra 
sources of uncertainties (e.g. , those due to the lack of 
precise knowledge of SM parameter values) and relevant 
external information (e.g. , about the parameter range) 
are difficult to accommodate in this scheme. 

In this work we report results from a Bayesian explo- 
ration of the CMSSM parameter space, obtained through 
the use of Markov Chain Monte Carlo methods. In partic- 
ular we focus on the prospects for direct neutralino dark 
matter detection with the next generation of dark matter 
searches. We refer the reader to [8] for full details. The 
Bayesian approach has several technical and statistical ad- 
vantages over the more traditional fixed-grid scan tech- 
nique, the most important being perhaps the ability to 
incorporate all relevant sources of uncertainties, e.g. the 
residual uncertainty in the value of SM parameters. This 
means that the inferred high probability regions in terms 
of e.g. neutralino mass and scattering cross section take 
fully into account all sources of uncertainty relevant to the 
problem. For other works applying a similar approach to 
the CMSSM, see [9] (with some relevant differences) and 
more recently [10,11]. 

2. Parameter space, priors and data used 

We restrict our analysis to the case sign(^) = +1, as 
motivated by the fact that the observed anomalous mag- 
netic moment of the muon is positive, and since the sign of 
the SUSY contribution to it is the same as the sign of \x. 
We thus consider the 8 dimensional parameter space (8, ip), 
where 8 is a vector of CMSSM parameters, 



8 = (m ,TOi/ 2 , A ,tan/3) 
while tp is a vector of relevant SM parameters, 



(1) 



iP = (M t ,m b (m b ) Mb ,a em (M z ) Mb ,a s (M z ) Mb ), (2) 

where M t is the pole top quark mass, mb(mb) MS is the bot- 
tom quark mass at nib, while a em (Mz) MS and a s (Mz) MS 
are the electromagnetic and the strong coupling constants 
at the Z pole mass Mz, the last three evaluated in the 
MS scheme. We are not interested in constraining the 
value of the SM parameters, but rather in including the 
effect of the uncertainty in their experimental determina- 
tion on the mass spectra, cosmological dark matter abun- 
dance and other observable quantities. This in turns has a 
non-negligible impact when making inferences on the high 
probability regions for the CMSSM parameters, 8, or for 
any other derived quantity. At the end of the analysis, the 



set of so called "nuisance parameters" ip is integrated out 
from our probability distribution function (pdf). It turns 
out that including the nuisance parameters in our analysis 
as free parameters (rather than fixing them to the central 
experimental value as it is done in most grid scans) has 
an important impact in widening the constraints on the 
CMSSM parameters, an effect that should not be ignored 
when carrying out proper statistical inference on 8 or on 
other quantities of interest such as the dark matter scat- 
tering cross section. 

Bayesian statistics makes use of Bayes theorem, 



P(v\d) 



p(d) 



(3) 



where we have introduced n = (8, tp), to compute the pos- 
terior probability distribution p(n | d). On the rhs ofEq. (3), 
the quantity p(d | n, f(n)), is called the likelihood and it sup- 
plies the information provided by the data, by comparing 
the base parameters n or any derived function f(n) to the 
data d. The quantity n(n) denotes a prior probability den- 
sity function (hereafter called simply a prior) which encodes 
our state of knowledge about the values of the parameters 
before we see the data. We take the prior to be flat (i.e., 
constant) in the variables n, and we further need to specify 
the range along each direction. If the constraining power 
of the likelihood is strong enough to override the choice of 
the prior, than the latter does not matter in the final infer- 
ence based on the posterior pdf. We found this to be the 
case for all of the parameters but Too and sparticles masses 
that mostly depend on it, such as sleptons and squarks. For 
those quantities, choosing a prior 50 GeV < m < 2 TeV 
cuts away a large region of parameter space that is not dis- 
favoured by data, the so called "focus point region" . If the 
prior range is enlarged to 50 GeV < mo < 4 TeV, then 
a considerable part of the focus point region is included 
in the analysis, even though it is still not possible to con- 
strain the value of mo from above. With presently available 
data, any upper limit on m is purely a consequence of the 
prior adopted. In the following we present results for the ex- 
tended range with a prior region up to 4 TeV for too, mi/2, 
and in the range \A Q \ < 7 TeV, 2 < tan/3 < 62. The prior 
range on the nuisance parameters does not influence the 
final results, since the SM parameters are rather tightly 
constrained by the data (summarized in the top section of 
Table I). 

For our analysis, from the CMSSM and SM parame- 
ters n we compute a series of derived observable quantities 
f(n): the W gauge boson mass, the effective leptonic weak 
mixing angle sin 2 8 eS , the anomalous magnetic moment of 
the muon, = (g — 2) M , the branching ratios BR(B — ► 
X s -f) and BR(B S — > /i+^~), the cosmological neutralino 
relic abundance f^cDM^ 2 (by solving the Boltzmann equa- 
tion numerically as in [12]), the light Higgs mass and the 
superpartner masses (computed with the package SOFT- 
SUSY vl.9 [13]). For all of those quantities, relevant mea- 
surements (summarized in the bottom section of Table 1) 
or experimental limits arc included via the likelihood and 



Nuisance parameter 


Mean value 


Uncertainty 






experimental theoretical 


M t 


172.7 GeV 


2.9 GeV 


N/A 




4.24 GeV 


0.11 GeV 


N/A 


a s (M z )MS 


0.1186 


0.002 


N/A 


l/aem(M Z )^ 


127.958 


0.048 


N/A 


Derived observable 




M w 


80.425 GeV 34 MeV 


13 MeV 


sin 2 e cB 


0.23150 


16 x 10" 5 


25 x 10" s 


fegUSY x 10 i() 


25.2 


9.2 


1 


BR(B -» X s7 ) x 10 4 


3.39 


0.30 


0.30 


fi x h 2 


0.119 


0.009 


o.i n x h 2 



Table 1 

Top section: experimental mean fi and standard deviation for the 
nuisance parameters used in the analysis. Bottom section: as above, 
but for derived observable quantities, including a theoretical uncer- 
tainty describing the imprecise mapping of CMSSM and SM param- 
eters onto observable quantities. Notice that the value for the cos- 
mological neutralino abundace il x h? stems from WMAP 1st year 
data combined with other cosmological observations. 

used to constrain high posterior probability regions of the 
model. We assume that all of the cold dark matter is in the 
form of ncutralinos. Furthermore, the likelihood is modi- 
fied in such a way that it includes estimated theoretical 
uncertainties in the mapping from CMSSM and SM pa- 
rameters to derived quantities, another major advantage 
of employing a Bayesian approach (see [8] for details). We 
then compute a spin-independent dark matter WIMP elas- 
tic scattering cross section on a free proton, a^ 1 , including 
full supersymmetric contributions which have been derived 
by several groups [14,15,16,17,18], but we do not include 
current constraints in the likelihood, in view of the uncer- 
tainties in the structure of the Galactic halo (e.g. , exis- 
tence of clumps of dark matter and therefore the value of 
the local halo mass density) as well as in the values of some 
hadronic matrix elements entering the computation of cr^ 1 . 
Our Bayesian approach allows to easily compute the poste- 
rior pdf for the cross section or any other derived variable. 
As we now discuss, our result for high-probability regions 
is highly suggestive of a possible detection by the next gen- 
eration of direct dark matter searches. 

3. Results for neutralino direct detection 

In Fig. 1 we present the 2-dimensional posterior pdf 
for (Tp 7 and m x , with all other parameters marginalized 
over. For comparison, we also show current CDMS-II [19], 
Edelweiss-I [20] and UKDMC ZEPLIN-I [21] 90% CL up- 
per limits, but we stress that this constraint has not been 
used in our analysis. 

In Fig. 1 the biggest, banana-shaped region of high prob- 
ability (68% regions delimited by the internal solid, blue 
curve) shows a well-defined anticorrelation between a^ 1 
and m x . It results from two allowed regions in the CMSSM 
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Fig. 1. The 2— dimensional probability density in the neutralino mass 
and spin— independent cross section plane in the CMSSM (all other 
parameters marginalized) with the contours containing 68% and 95% 
probability also marked. Current 90% experimental upper limits 
are also shown. A large fraction of the high— probability region lies 
just below current constraints and it will be probed by the next 
generation of dark matter searches, starting from the focus point 
region (horizontal region at of 7 ~ 10~ 8 ). 



parameter space: the bulk and stau coannihilation region 
and from the A-resonance at large tan/3. This region cov- 
ers roughly the range 10~ 10 < 1 < 10~ 8 pb and 200 < 
m x < 700 GeV. In both cases the dominant contribution 
to dp 1 comes from a heavy Higgs exchange. At small m x < 
100 GeV we notice a small vertical band of fairly low proba- 
bility density (< 0.2) at small a^ 1 , a region where the light 
Higgs resonance contribute to reducing fl x h 2 at small miy 2 
to meet the WMAP measurement. It is interesting to no- 
tice that our Monte Carlo procedure is able to resolve such 
small features of parameter space, that are almost invari- 
ably missed by usual fixed-grid scans. This region would 
also disappear with a fair improvement in the lower bound 
on nih ■ Finally, we can see a well pronounced region of high 
probability at fairly constant 1 ~ 1.6 x 10~ 8 pb for m x < 
420 GeV which at low m x partly overlaps with the previous 
region. At 95% this region extends up to m x < 720 GeV for 
fairly constant a^ 1 . This "high" a^ 1 band results from the 
focus point region, basically independently of tan/3. This 
result has to be interpreted carefully, since there are large 
uncertainties associated with FP region, in particular with 
its location in the (m 1 / 2 ,rn n ) plane mentioned earlier. De- 
spite those outstanding questions, we believe that it is safe 
to expect that the FP will be the first to be probed by dark 
matter search experiments. 

After marginalizing over all other parameters, we obtain 
the following 1-dimensional regions encompassing 68% and 
95% of the total probability (see Fig. 2): 
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Fig. 2. The 1-dimensional probability distribution for the 
spin— independent neutralino scattering cross section in the CMSSM. 
The three curves explore the impact of different assumptions in the 
analysis: assuming a 4 TeV upper prior on mo (solid), a 2 TeV prior 
(dotted) or removing from the analysis the constrain from g — 2 
(dashed, for the 4 TeV prior case) The vertical, dashed lines delimit 
the region encompassing 95% probability for the 4 TeV prior and 
including g — 2. The spike around a^I ^ 10 -8 corresponds to the fo- 
cus point region at large mo. Future one— tonne detectors will reach 
down to <jSI ^ 10~ 10 and thus explore most of the 95% region. 

1.0 x 10~ 10 pb < cr^ < 1.0 x 10~ 8 pb (68% region), 
0.5 x 10~ 10 pb < cr| J < 3.2 x 10~ 8 pb (95% region). 

(4) 

Currently running experiments (most notably CDMS-II 
but also Edelweiss-II and ZEPLIN-II) should be able to 
reach down to a few x 10~ 8 pb, on the edge of exploring this 
FP region. A future generation of "one-tonne" detectors 
is going to reach down to a^ 1 > 10 _10 pb, thus exploring 
almost the whole 68% region and much of the 95% interval 
as well. 

This result for the pdf on is rather robust with re- 
spect to a change in the prior range. The spike at a^ 1 ~ 
10~ 8 disappears if one cuts away most of the focus point by 
imposing a prior tuq < 2 TeV, but the 95% region is only 
slightly shifted to lower values. We have also considered the 
impact of removing from the analysis the constrain coming 
from g — 2. We have found that in this case mo becomes es- 
sentially unconstrained, and this adds statistical weight to 
the focus point region. However, even in this case the high- 
probability range for a~ remains close to the values given 
in Eq. (4), shifting to 0.3 x 10~ 10 pb < trf 1 < 2 .8 x 10" 8 pb 
(95% probability region). 

We notice that the high-probability regions as described 
above do not necessarily coincide with the best fitting 
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Fig. 3. Distribution of the values of a^ 1 in the (m]/ 2 ,mo) plane. The 
dots delimit the 95% confidence region in parameter space. Values 
down to o-p 1 ~ 10 — 8 (darker dots) will be among the first to be 
explored by currently running dark matter searches. This includes a 
large parte of the focus point region at large mo. 



points in parameter space if the pdf is strongly non- 
Gaussian, as in the present case. We refer the reader to [8] 
for a detailed description of the discrepancy and a discus- 
sion of its meaning in terms of probabilistic inference. 

Finally, it is interesting to consider the complementar- 
ity between direct dark matter searches and the reach of 
future collider experiments, especially in view of the immi- 
nent start of operations at LHC. We plot in Fig. 3 the distri- 
bution of values of a^ 1 in the plane (rrii/2, mo), where the 
points roughly cover the region encompassing 95% of prob- 
ability. It is apparent that the largest values of the cross 
section (<7p > 10~ 8 ) are clustered in the region 200 GeV < 
m i/2 ^ 500 GeV, covering values of mo all the way up to 
our prior limit of 4 TeV in the focus point region. As the 
currently running dark matter detectors start probing value 
down to <7p 7 ~ 10 -8 , they will begin exploring the focus 
point region of parameter space, which lies at large values 
of mo that are largely outside the reach of LHC searches. 
From our analysis a very promising synergy thus emerges 
between direct dark matter collider searches, that will soon 
begin to squeeze from different directions the current high- 
probability region of the CMSSM parameters space. 

4. Conclusions 

We have presented a detailed investigation of the 
prospects for dark matter detection in the framework of 
the CMSSM parameter using state-of-the-art Bayesian 
methods. The power and flexibility of the approach allows 
to probe many previously unexplored ranges of parameters 
and to fully incorporate the effects of remaining uncer- 
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tainties in relevant SM parameters and other theoretical 
uncertainties in computing mass spectra and observables. 

We have shown that the WIMP dark matter direct de- 
tection elastic scattering cross section 1 presents a wide 
spread of values (below today's limits) at around 10 _9±1 pb 
and a strong anticorrelation with m x . In addition, a region 
at relatively large a^ 1 ~ 1.6 x 10~ 8 pb, and fairly inde- 
pendent of m x , appears to be a feature of the focus point 
region (despite large theoretical uncertainties) and will be 
the first to be tested in direct detection experiments. 

Despite the remaining uncertainties arising from the lo- 
cal dark matter density profile and the large errors on nu- 
clear form factors, it seems that a large fraction of the 
high-probability parameter space of neutralino dark mat- 
ter in the CMSSM will be within reach of currently running 
upgraded dark matter detector, and most of it will be ex- 
plored by future one-tonne detectors. We have highlighted 
the complementarity of these searches with direct SUSY 
searches at colliders and in particular showed that the fo- 
cus point region, mostly outside the reach of the LHC, will 
be among the first to be probed by direct dark matter de- 
tection experiments. Those results are largely robust with 
respect to changes in the a priori allowed range for the 
parameters or to the exclusion of the anomalous magnetic 
moment of the muon measurement from the analysis. 
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